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Resumen 


El presente estudio analizó la influencia del cambio climático en el comportamiento 
de los caudales punta en la Región Metropolitana, Región del Libertador Bernardo 
O'Higgins y la Región del Maule, en la zona central de Chile. Para ello, se utilizó 
la FDP de Gumbel, con la que se obtuvieron los mejores resultados para ajustar los 
caudales punta, los cuales se analizaron para los periodos de retorno de 10, 20, 30, 40 
y 50 años. Se concluye, gracias a la segmentación de los datos a través de series en un 
periodo dado, que los caudales presentan una tendencia a registrar mayores valores 
en los periodos más actuales. Esto se visualizó con los valores promedio, gráficas 
de tendencia y ajustes de funciones de distribución de probabilidad. Finalmente se 
pudo evidenciar de manera gráfica y estadística que en los últimos años ha existido 
un aumento de los caudales máximos anuales, lo cual podría confirmar desde el 
punto de vista hidrológico procesos de cambio climático que de forma eventual han 
generado un mayor aporte de caudales a través de la disminución de las reservas de 
nieve y el retroceso glaciar por el aumento de las temperaturas. 


Palabras clave: caudales punta, cambio climático, funciones de distribución de 


probabilidad, zona mediterránea. 


Introducción 


Existen evidencias de que en los últimos años 
el planeta ha experimentado alzas en sus 
temperaturas, y debido a que Chile se encuentra 
en una zona de transición climática, el fenómeno 
del cambio climático podría ser muy impactante, 
en particular para la disponibilidad y el uso 
del recurso hídrico. Por ello, la influencia de 
este fenómeno sobre el comportamiento de los 
caudales máximos podría verse aumentada 
(Ivanova y Corredor, 2006; Fuenzalida et al., 
1989). 

Un estudio realizado en California, Estados 
Unidos, asociado con el impacto del cambio 
climático, enfocado desde el punto de vista 
de la hidrología, consistió en el análisis y la 
aplicación de modelos para determinar la 


sensibilidad de algunas variables hidrológicas, 
donde se encontraron grandes cambios en la 
magnitud de los caudales punta, obteniéndose 
un mayor aumento de éstos en proporción a las 
precipitaciones; además, se pudo demostrar en 
este estudio, que el aumento de las temperaturas 
estarían derritiendo o disminuyendo los sectores 
con presencia de nieve (Miller et al., 2003), misma 
situación que podría estar ocurriendo en la zona 
mediterránea de Chile. 

Las variables hidrológicas, como las preci- 
pitaciones y los caudales, muestran un papel 
fundamental en cuanto a su sensibilidad al 
cambio climático, generando inquietudes sobre 
su comportamiento, debido a que son variables 
directamente influenciables por fenómenos de 
amplio espectro temporal y espacial (Whitaker 
et al., 2002; Waylen y Woo, 2008). 


Tecnología y Ciencias del Agua, vol. IV, núm. 2, abril-junio de 2013, pp. 5-19 


Ciencias del Agua, vol. IV, núm. 2, abril-junio de 2013 


«o180/0u3aL 


Pizarro et al., Influencia del cambio climático en el comportamiento de los caudales máximos en la zona Mediterránea de Chile 


Las crecidas son procesos naturales, sin 
periodicidad, constituidos por un incremento 
importante y repentino del caudal en un 
sistema fluvial, el cual lleva consigo un ascenso 
del nivel de la corriente, que puede desbordar 
el cauce menor para ocupar progresivamente 
el cauce mayor hasta alcanzar un máximo 
o caudal-punta y descender a continuación 
(Aparicio, 1997; Ollero, 1996; Chow et al., 1994). 
Según Paoli et al. (1998), las crecidas que se 
presentan en términos hidrológicos poseen 
un grado de riesgo y una probabilidad de 
excedencia diferente, según el caudal máximo, 
el volumen o la duración que se considere. 

Dicho lo anterior, y con el fin de entender 
las posibles influencias del cambio climático 
sobre los caudales punta de la zona central de 
Chile, es necesario analizar el comportamiento 
de los caudales máximos debido a que ellos son 
un reflejo de las situaciones particulares que 
pueden estar sucediendo en la cuenca (Heusser 
et al., 1981). 

El aumento de las frecuencias de inun- 
daciones, la escasez de agua para millones 
de personas, la amenaza del incremento 
del nivel del mar y la destrucción de 
comunidades costeras y hábitats (Combes 
et al., 2003; Kroll y Vogel, 2002; Sullivan et 
al., 2007) son preocupaciones que el cambio 
climático puede traer consigo. En este sentido, 
estudios realizados por Carrasco et al. (2005 
y 2008), analizaron el comportamiento de 
la isoterma 0”, variable climática que afecta 
de manera 
nivoglaciar, encontrándose que dicha variable 


directa al comportamiento 
ha aumentado en su altura, es decir, se ha 
producido un aumento de la temperatura 
(Rosenbluth et al., 1997), y por consiguiente, 
el volumen que anteriormente era nieve, a 
raíz del cambio de cota de la isoterma, se 
transforma de forma directa en escorrentía 
superficial. 

Por tanto, surge la necesidad de determinar 
si los caudales máximos de cuencas de la 
zona central de Chile han presentado alguna 
variación significativa en las últimas décadas, 
producto de un aparente cambio climático. 


Antecedentes generales 


El estudio abarcó las regiones Metropolitana, 
del Libertador Bernardo O'Higgins y del 
Maule, ubicadas en la zona central de Chile, 
caracterizadas por un clima semiárido oO 
subhúmedo, con precipitaciones que van desde 
los 250 mm hasta los 500 mm anuales, y con 
sistemas hidrográficos de régimen pluvial y 
pluvionival. 

La Región Metropolitana se ubica entre 
los 32% 557 y los 34” 19 de latitud sur, y 
entre los 69” 46' y 71” 43” de longitud oeste, 
aproximadamente. Esta región limita al sur 
con la Región del Libertador General Bernardo 
O'Higgins. Tiene una superficie de 15 403.2 
km”, lo que representa un 2.05% del territorio 
nacional (excluido el Territorio Chileno 
Antártico). 

Por otro lado, la Región del Libertador 
Bernardo O'Higgins está situada geográ- 
ficamente entre los paralelos 33* 51' y 35” 01' de 
latitud sur, y 72” 02” y 70” 02” de longitud oeste. 
La superficie territorial es de 16 387 km, lo 
que en relación al total de Chile continental en 
superficie corresponde al 2.3%. Limita al norte 
con la Región Metropolitana. 

La Región del Maule se extiende entre los 34” 
41” y 36” 33” de latitud sur, y desde los 70” 20” 
de longitud oeste hasta el océano Pacífico. La 
superficie calculada por el Instituto Geográfico 
Militar es de 30 296 km?, equivalente al 4% 
del territorio nacional. Al norte limita con la 
Región de O”Higgins, al sur con la Región del 
Bío-Bío, al oeste con el océano Pacífico y al este 
con la cordillera de los Andes. 

En la figura 1 se presenta la ubicación 
de las estaciones fluviométricas por región. 
Asimismo, en el cuadro 1 se tienen los nombres 
de las estaciones, periodo de datos de cada una 
de ellas, caudal promedio y años de datos. 


Metodología 
La información necesaria para la realización 


del presente estudio fue proporcionada por 
la Dirección General de Aguas (DGA), y ésta 
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Figura 1. Regiones Metropolitana de O'Higgins y del Maule, y sus estaciones fluviométricas. 


correspondió a los caudales máximos anuales 
de seis estaciones fluviométricas de la Región 
Metropolitana, 16 estaciones de la Región de 
O'Higgins y 28 estaciones de la Región del 
Maule (cuadro 1). 

La información recopilada fue entregada 
por el organismo hidromensor oficial de Chile, 
la Dirección General de Aguas, por lo que no 
fue necesario realizar transformaciones con 
limnigramas, ni curvas de descarga, para 
obtener los caudales máximos anuales. Luego, 
para cada serie de datos, se determinaron, por 


una parte, los estadígrafos de posición, también 
llamados de tendencia central, para indicar 
alrededor de qué valor se agruparon los datos 
obtenidos, como es el caso de la media, y por 
otra parte, los estadígrafos de dispersión. Estos 
estadígrafos o estadísticos extraen información 
de una muestra, indicando las características 
de la población. En el caso del presente estudio, 
los estadígrafos utilizados fueron la media, 
cuyo cálculo es el valor esperado de la variable 
misma O primer momento respecto al origen 
(muestra la tendencia central de la distribu- 
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Cuadro 1. Nombre de las estaciones fluviométricas, periodo de registro de datos, años de la serie y caudal promedio 


de cada estación. 
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Región Núm. Estación Periodo Jos ta a 
E 1 Río Maipo en El Manzano 1979-2007 43 439.78 
E 2 Río Maipo en El Cabimbao 1965-2007 20) 1 443.5 
al 3 3 Río Maipo en Las Melosas 1965-1972*,1974-2007 27 162.42 
2 Z 4 Río Maipo en San Alfonso 1965-2007 43 267.45 
El 5 Río Mapocho en Los Almendros 1965-1968*, 1979-2007 33 84.41 
e 6 Río Mapocho en Rinconada de Maipú 1965-2007 43 253.57 
yd Canal Sauzal en Puente Termas 2005-2007 59.47 
8 Canal Teno 2002-2003 39.9 
9 Estero Alhué en Quilamuta 1970-2007 38 172.51 
10 Estero Nilahue en Santa Teresa 1985-2007 23 456.06 
n Estero Zamorano en Puente El Niche 1985-2007 23 543.49 
el 12 Río Cachapoal 5 km aguas abajo Junta Cortaderal 1989-2007 19 486.8 
-] 13 Río Cachapoal en Puente Arqueado 2002-2007 5 113.51 
o 14 Río Cachapoal en Puente Termas de Cauquenes 2002-2007 9 68.96 
3 15 Río Claro en El Valle 1970-2007 38 223.03 
S 16 Río Claro en Hacienda Las Nieves 1960-2007 48 109.08 
Z 17 Río Claro en Tunca 1994-2007 14 775.21 
18 Río Cortaderal ante Junta Río Cachapoal 1985-1991*/2006-2007 24.96 
19 Río Las Leñas ante Junta Río Claro 2006-2007 11.13 
20 Río Pangal en Pangal 1985-2007 23 99.63 
21 Río Tinguiririca Bajo Los Briones 1984-2007 23 378.61 
22 Río Tinguiririca en Los Olmos 2002-2007 6 79.28 
23 Río Achibueno en La Recova 1987-2006 20 46.19 
24 Río Ancoa en El Morro 1960-2006 47 22,13 
25 Río Loncomilla en Las Brisas 1983-2005 23 257.49 
26 Río Claro en Camarico 1963-2006 42 Tata 
27 Río Cauquenes en Desembocadura 1986-2006 21 15.12 
28 Río Colorado en Junta con Palos 1975-2006 82 45.78 
29 Río Claro en Los Queñes 1986-2006 21 20.56 
30 Río Maule en Armerillo 1960-1978*, 2001-2006 20 198 
31 Río Cauquenes en El Arrayán 1987-2006 20 8.78 
32 Río Loncomilla en Bodega 1985-2006 22 145.16 
33 Río Longaví en El Castillo 1964-2006 42 34.73 
E] 34 Río Maule en Forel 1988-2006 16 494.65 
E 35 Río Mataquito en Licantén 1987-2006 20 132.72 
a 36 Río Loanco en Desembocadura** 1987-2006 20 4.72 
E 37 Río Maule en Logitudinal 1962-2005 43 204.67 
bo 38 Río Melado en El Salto 2003-2006 8 102.81 
E 39 Río Purapel en Nirivilo 1960-2006 47 2.73 
40 Río Perquilauquén en Ñiquen 1987-2006 20 39.66 
41 Río Palos en Junta con Colorado 1967-2006 40 29.14 
42 Río Putagán en Yerbas Buenas 1987-2006 20 21.62 
43 Río Perquilauquén en Quella 1963-2005 43 57.74 
44 Río Longaví en Quiriquina 1960-2006 47 46.62 
45 Río Lircay en Puente Las Rastras 1961-2006 45 13.88 
46 Río Claro en Rauquén 1999-2006 8 95.83 
47 Río Perquilauquén en San Manuel 1960-2006 47 33.05 
48 Río Teno después de Junta con Claro 1960-2006 47 58.19 
49 Río Teno Bajo Quebrada Infiernillo 1985-2006 19 32.62 
50 Río Purapel en Sauzal 1999-2005 Y) 4.97 


* Estaciones con información intermedia faltante; ** estación perteneciente a las cuencas costeras del Maule, en el límite de la región. 
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ción); la varianza, que mide la variabilidad de 
los datos (es el segundo momento, con respecto 
al origen), y finalmente la desviación estándar, 
que es una medida de la variabilidad, ya que 
es la raíz cuadrada de los cuadrados de las 
diferencias. 

Posteriormente, los datos fueron ajustados 
a cuatro funciones de distribución de proba- 
bilidad (FDP), a saber: Función de Distribución 
Log-Normal, Pearson Tipo Ill, Goodrich y 
Gumbel, con el fin de determinar cuál de estas 
funciones se adecuaba más a los caudales 
máximos anuales. Dichas funciones se definen 
en el cuadro 2. 

Es imprescindible considerar que la infor- 
mación de parámetros estadísticos se calculó 
con la totalidad de datos proporcionados, a 
fin de no realizar ningún tipo de completación 
de datos faltantes, lo que le da un carácter 
fidedigno a la información utilizada. En los 
casos en que faltó información para uno o más 
años intermedios, se optó por no completar, 
para no producir sesgos de ninguna forma en 
la información original. 

Después se determinó la calidad de las 
FDP ajustadas, en donde se realizó una 
contrastación para cada dato que se obtuvo en la 
frecuencia acumulada y en la frecuencia teórica 
acumulada; para ello se utilizaron métodos 
cuantitativos, a fin de determinar la bondad del 
ajuste, como el coeficiente de determinación R? 
y el test de Kolmogorov-Smirnov. 

Para conocer el valor a asumir por la 
variable aleatoria de caudales máximos, se 
utilizaron los periodos de retorno de 10, 20, 30, 
40, 50 y 100 años. 


Análisis y discusión de resultados 


A continuación se presentan los principales 
resultados obtenidos en el presente estudio, 
acompañados de su discusión y análisis. 


Calidad de los ajustes 


De acuerdo con los resultados alcanzados de 
bondad de ajuste y por la facilidad del cálculo, 
se concluyó que en las tres regiones, la función 
de distribución de probabilidad de Gumbel 
es la más adecuada para su aplicación a los 
caudales máximos (Yue et al., 2006). 

A continuación se presentan los resultados 
entregados por ambas pruebas en cuanto a 
la calidad del ajuste de las series de caudales 
máximos de las tres regiones en estudio 
(cuadro 3). 

Por lo tanto, se puede determinar que las 
FDP Gumbel y Log-Normal, presentan los 
valores más altos de R? (cercanos al 95%), y 
un 100% de aceptación de K-S, por lo que se 
determinó que las dos funciones son adecuadas 
para el ajuste de los caudales máximos anuales. 
Sin embargo, y aunque la Distribución Log 
Normal presenta un R* promedio levemente 
más alto que el resto de las funciones, se utilizó 
la distribución de Gumbel debido a que es la 
más fácil de determinar y a la vez presenta una 
menor cantidad de parámetros a estimar. 


Variabilidad regional 


Con el fin de comprender la relación y el 
comportamiento existente entre los caudales 


Cuadro 2. Funciones de distribución de probabilidad utilizadas. 


FDP Log-normal FDP Pearson Tipo IM 


FDP de Goodrich FDP de Gumbel 


F(x) = 


-0 


_1/Inx-a 
2 


di las 


La j Al 
El) =_——— fe ar 
(a) al (B) + 
x: representa el valor a asumir la variable aledtória 
aL, P y Ó: parámetros. 


e: constante de Neper. 


por la variable aleatoria. 
P, x: parámetros. 
e: constante de Neper. 


l) 


x: representa el valor a asumir por 


F(X)=1- pa) 


d(x-4) 


FEO) =e-e 
x: representa el valor a 
asumir por la variable 
aleatoria. 


x: representa el valor a 
asumir por la variable 

aleatoria. 

e: constante de Neper. 

u y d: parámetros. 


a, p y X,: parámetros. 
e: constante de Neper. 
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Cuadro 3. Bondad promedio del R? y porcentaje de aceptación promedio K-S para el conjunto de estaciones, 


para las tres regiones en estudio. 


Función de Distribución de Probabilidad (FDP) 
Gumbel Log-Normal Goodrich Pearson tipo III 
Promedio R?(%) 95.45% 98% 69.30% 95.40% 
Región 
Metropolitana a de 100% 100% 0% 100% 
aceptación K-S 
Promedio R?(%) 95.45% 93.88% 87.75% 94.92% 
Región de O'Higgins : 
P taje d 
o 100% 100% 100% 100% 
aceptación K-S 
Promedio R?*(%) 96.40% 94.10% 97.20% 94.90% 
Región del Maule ¡ 
: Earenajede 100% 100% 96.4% 82.1% 
aceptación K-S 


que se presentaron en este estudio, se elaboró 
un índice relativo de magnitud de crecidas 
para las tres regiones, el cual representa la 
relación existente entre la estación que presenta 
el mayor caudal y las demás estaciones con 
distintos valores de caudal. 

En este índice se consideraron los valores 
reales empíricos y los valores para periodos de 
retorno T = 30 y T'= 50, el cual queda definido 
con la siguiente expresión: 


mes QmáxEmayor 


Omáxi 


Donde: 

IRmc = índice relativo de magnitud de crecidas. 
Qmáx Emayor = caudal máximo anual de la 
estación con el valor mayor. 

Omáx ¡ = caudal máximo, estación 1. 

Para la Región Metropolitana se consideró 
como estación patrón a la estación El Cabimbao; 
para la Región del Libertador Bernardo 
O'Higgins lo fue Santa Teresa, y para la Región 
del Maule lo fue Río Maule en Forel. 

El cuadro 4 y la figura 2 presentan los 
resultados del índice de magnitud de crecidas 
tanto para los valores reales de caudales 
máximos como para los valores probables en 
los periodos de retorno de 30 y 50 años, con el 


objetivo de entender sus comportamientos a 
través del tiempo. 

Como se puede observar en el cuadro 4 y la 
figura 2, la estación Río Maipo en El Cabimbao 
sobrepasa 15.4 veces el caudal máximo real de 
la estación Río Mapocho en Los Almendros, 
donde esta estación es la que presenta los 
valores más bajos y se mantiene a través del 
tiempo tanto para 30 como para 50 años. 

En cuanto a la Región de O'Higgins, los 
valores observados muestran que la estación 
Río Nilahue en Santa Teresa es la que presenta 
un mayor valor de caudal máximo, y al ser com- 
parado con la estación Río Pangal en Pangal, 
ésta la supera en 7.4 veces, disminuyendo en T' 
=50 a 6.9 veces. 

Para el caso de la Región del Maule, la 
estación Río Maule en Forel sobrepasa 54.6 
veces el caudal máximo de la estación Río 
Purapel en Nirivilo, siendo esta última la que 
presenta los valores de caudales máximos 
más reducidos entre las estaciones, lo que 
se mantiene en el tiempo tanto para 30 como 
50 años, llegando a sobrepasar 69.1 veces los 
caudales máximos probables en un periodo de 
retorno de 50 años. Sin embargo, la estación Río 
Loncomilla en Las Brisas sólo es sobrepasada 
2.3 veces por la estación Río Maule en Forel. 

En resumen, se puede observar que el 
caudal máximo en las tres regiones en estudio 
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Cuadro 4. Caudales máximos reales para estaciones de la Región Metropolitana, la Región de O'Higgins 


y la Región del Maule. 
Caudal 009 codo di 
Estaciones Máximo no oa 
(m/s) tio pa e T=30 | T=50 
(T = 30) (T= 50) di 
po Río Maipo en El Cabimbao 5 095.4 4 387.1 4 929.9 il il 1 
E Río Maipo en El Manzano 1 430.3 1 044.3 1 155.7 3.6 4.2 4.3 
S Río Mapocho en Rinconada de Maipú 1 124.7 780.6 877.8 4.5 5.6 5.6 
S Río Maipo en San Alfonso 740.4 589.2 648.5 6.9 7.4 7.6 
Sl Río Maipo en Las Melosas 427.7 340.6 373.4 11.9 1218) 13.2 
Ez Río Mapocho en Los Almendros 330.2 294.6 333.3 15.4 14.9 14.8 
Estero Nilahue en Santa Teresa 1 396.7 1298,9 1 454.3 ll dl 1 
Estero Zamorano en Puente El Niche 1 189.1 1 480.1 1 652.8 1,2 0.9 0.9 
a Río Tinguiririca Bajo Los Briones 1 007.3 945.6 1 050.1 1.4 1.4 1.4 
E Río Claro en Tunca 775.2 879.4 980.6 1.8 15 15 
e Río Claro en El Valle 621.1 564.6 627.5 22) za ZiS 
E Estero Alhué en Quilamuta 547.2 540.7 608.6 2.6 2.4 2.4 
Es Río Cachapoal 5 km aguas abajo Junta Cortaderal 486.8 436 475.6 E) 3 3.1 
Río Claro en Hacienda Las Nieves 451.4 345.2 388.7 3.1 3.8 3.7 
Río Pangal en Pangal 188.5 192.2 209.3 7.4 6.8 6.9 
Río Maule en Forel 17 212.9 17 734.6 20 039 10 10 10 
Río Loncomilla en Las Brisas 7 623.4 77434 8 498.6 DEa) 2.3 24 
Río Mataquito en Licantén 4 638.3 4 844.1 5 324.4 3,7 37 3.8 
Río Loncomilla en Bodega 4 227.8 4 524.3 5028 4.1 3.9 4.0 
Río Maule en Logitudinal 3064 3061.1 3 371.3 5.6 5.8 5/9 
Río Longaví en El Castillo 2877.9 2900 3 184.6 0.9) 6.1 6.3 
Río Achibueno en La Recova 2 621.2 2 552.6 2750 6.6 69 18 
Río Longaví en Quiriquina 2 534.9 2 482 2 740.8 6.8 71 7.3 
Río Perquilauquén en Quella 2 319.6 2 420.7 2 665.3 74 7.3 75 
Río Maule en Armerillo 2 305.2 2 398 2 637.7 75 7.4 7.6 
Río Claro en Rauquén 2 100.1 2 309.7 2 560.8 8.2 7.7 7.8 
3 Río Claro en Camarico 1 723.8 1593.5 1745.5 no) 11.1 ES 
S Río Perquilauquén en Ñiquen 1 378.6 1196.6 1317.4 125 14.8 15,2 
E Río Perquilauquén en San Manuel 1 244.9 1184.5 1 303.5 13.8 15.0 15.4 
sl Río Teno después de Junta con Claro 1 155.8 1 083.6 1 205.6 14.9 16.4 16.6 
3) 
pá Río Ancoa en El Morro 1 080.9 1 033.3 1 139.6 15.9 17.2 17.6 
Río Cauquenes en El Arrayán 1 032.1 1 025.4 1 126.6 16.7 17.3 17.8 
Río Melado en El Salto 972 1 010.2 1 042.2 17.7 17.6 19.2 
Río Cauquenes en Desembocadura 922.4 915.5 1 029.1 18.7 19.4 19.5 
Río Loanco en Desembocadura** 854.9 780.2 857.6 20.1 227 23.4 pe] 
Río Lircay en Puente Las Rastras 850 712 787.5 20.3 24.9 25.4 S 
Río Colorado en Junta con Palos 689.6 658.1 ZA 24.9 26.9 27.4 ES 
Ria Palos en Junta con Colorado 597.4 563.9 619.5 28.8 31.4 32.3 A 
Río Claro en Los Queñes 527.9 558.7 616 32.6 31.7 32.5 Z 
Río Teno Bajo Quebrada Infiernillo 358.8 447.1 497.1 47.9 39.7 40.3 E 
Río Purapel en Sauzal 321.2 397.3 391.4 53.6 49.6 p1.2 £ 
Río Purapel en Nirivilo 315.2 261.2 289.9 54.6 67.9 69.1 > 
Donde: E 
Caudal máx: son los caudales reales máximos de cada serie. 3 
Caudal máx (T = 30): son los caudales máximos probables para un periodo de retorno de 30 años. - 
Caudal máx (T = 50): son los caudales máximos probables para un periodo de retorno de 50 años. = 
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Figura 2. IRmc caudales máximos anuales (T =30 y T =50 años), para las regiones Metropolitana, de O'Higgins y del Maule. 


posee un aumento gradual en el tiempo, cuya 
proyección ya da una idea de lo que está 
sucediendo con la variable en estudio. 


Años con más frecuencias en que existen 
caudales punta 


Con el fin de determinar si los caudales 
punta se han incrementado en las regiones en 
estudio en las últimas décadas, se confeccionó 
el cuadro 5, el cual presenta los periodos de 
tiempo y el número de caudales máximos que 
se han registrado en ese lapso en cada una de 
las estaciones de cada región. 

Para la Región Metropolitana se consideró 
el total de estaciones presentes (6), donde 
se muestra que los valores más altos fueron 
encontrados en la década de 1960. Para el 
caso de la Región del Libertador Bernardo 
O'Higgins, se contemplaron las nueve esta- 
ciones en estudio, de las cuales prácticamente 
la mitad presenta máximos entre los años 
2000-2006. En tanto, para la Región del Maule, 
se consideraron las 12 estaciones con mayor 
cantidad de registros, donde también se 


observa que la mayor cantidad de máximos se 
presenta en los últimos seis años. 

A fin de entender de mejor forma el com- 
portamiento de los caudales punta obser- 
vados anteriormente, se graficaron las tres 
estaciones de cada región que poseían la 
mayoría de los datos, con el fin de observar 
el comportamiento de los caudales punta, 
y adicionalmente las series de datos se divi- 
dieron en dos periodos. Así, la serie total, que 
corresponde a los años 1965-2006, se separó 
en dos series de datos; la primera consideró 
los años 1965-1981 y la segunda, 1982-2006. 
Por otra parte, la serie que considera el total 
de años desde 1970-2006, se separó en los años 
1970-1981 y 1982 2006. 

Se puede observar en un análisis más 
detallado de los caudales máximos en la 
estación Río Maipo en El Manzano, que el 
promedio general es de 448 m?/s, mientras que 
al separar las series, el primer periodo tiene un 
promedio de 351 m?/s, y en el caso del segundo 
periodo, alcanza un valor promedio de 513 
m/s. Esto queda demostrado en la línea de 
tendencia de los datos, que si bien no es tan 
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Cuadro 5. Número de estaciones por región con mayor frecuencia de existencia de caudales punta. 


De da Metropolitana B. O'Higgins Maule 
(núm. de estaciones) (núm. de estaciones) (núm. de estaciones) 
1960-69 3 - 
1970-79 1 1 2 
1980-89 - - 1 
1990-99 1 4 3 
2000-06 1 4 6 
Total estaciones analizadas 6 9 1122 


pronunciada, entrega valores más altos en los 
últimos años (figura 3). 

Para la estación Río Maipo en San Alfonso, 
se puede observar que los caudales máximos 
presentados para la primera estación presentan 
un promedio general de 283 m*/s. Al observar 
el primer periodo, el promedio alcanza 230 
m*/s, en tanto que el promedio del segundo 
periodo es levemente superior, con 308 m*/s. 
Esto queda demostrado en la línea de tendencia 
general, que también indica ascenso, aunque 
los datos son muy fluctuantes alrededor de la 
media. 

Finalmente, para la estación Río Mapocho 
en Hacienda de Maipú, el promedio general 
de la estación es de 281 m?*/s. Para el caso del 
promedio del primer periodo, sólo alcanza 
los 94 m?*/s, periodo que presenta una forma 
levemente ascendente. En el caso del segundo 


periodo, el promedio alcanza los 370 m*/s. Con 
esto se verifica que en los últimos años han 
aumentado los caudales, lo que claramente se 
evidencia con la línea de tendencia de toda la 
serie, la cual es muy pronunciada. 

En la estación Río Claro en Hacienda Las 
Nieves se presentan más fluctuaciones que en 
la estaciones anteriormente analizadas, donde 
el promedio general es de 120 m?/s. Para el 
primer periodo, el caudal promedio alcanza 
los 108 m?*/s; en el caso del segundo periodo, 
el promedio es levemente superior a la media y 
alcanza un valor de 129 m?*/s (figura 4). 

Para la estación Estero Alhué en Quilamuta 
se presentan muchas fluctuaciones con muchos 
valores extremos tanto bajos como altos; el 
promedio total alcanza los 177 m/s. En cuanto 
al primer periodo, el promedio denota un 
monto de 110 m?*/s, mientras que en el segundo 
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Figura 3. Caudales máximos, estaciones Río Maipo en El Manzano, Río Maipo en San Alfonso y Río Mapocho 


en Rinconada de Maipú, Región Metropolitana. 
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periodo, el promedio de los caudales máximos 
es de 208 m/s. 

Para el caso de la estación Río Claro en El 
Valle, el promedio general alcanza un total de 
250 m*/s, mientras que en el primer periodo 
el promedio es mayor, alcanzando los 228 
mí /s; este periodo cuenta con un par de datos 
extremos muy marcados. Para el segundo 
periodo, los valores de caudal son más o menos 
constantes, aumentando sólo en los últimos 
años; este periodo presenta un promedio de 
218 m?/s, valor menor a la media general. 

Para la estación Río Ancoa en El Morro, 
también sus valores son fluctuantes, pero se 
acercan más a la media, donde el promedio 
general alcanza los 317 m*/s. Asimismo, para 


el primer periodo de la serie, sigue siendo 
menor con 260 m?/s, y en el caso del segundo 
periodo, el caudal promedio es levemente 
superior a la media y alcanza un valor de 
356 m?*/s. Esta gráfica muestra una línea de 
tendencia general más pronunciada que las 
anteriores (figura 5). 

En el caso de la estación Río Purapel en 
Nirivilo, el promedio total es de 115 m* /s. En 
el primer periodo de datos considerado, este 
promedio es de 90 m?*/s, mientras que en el 
segundo periodo el valor asciende a 127 n?/s. 
Esta gráfica muestra una línea de tendencia 
general ascendente, pero no de forma 
pronunciada. Esta diferencia puede deberse al 


aumento progresivo de las temperaturas en las 
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Figura 4. Caudales máximos, estaciones Río Claro en Hacienda Las Nieves, Estero Alhué en Quilamuta 
y Río Claro en El Valle, Región del Libertador Bernardo O Higgins. 
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últimas décadas, tendencia que se observa en 
la gran mayoría de las estaciones analizadas. 

De las gráficas se puede inferir que los 
segundos periodos considerados muestran 
un valor promedio de caudal mayor que los 
primeros periodos; asimismo, las gráficas de 
las curvas de tendencia general son todas 
ascendentes, hecho que está claramente 
influenciado por los valores de caudales del 
segundo periodo. Además, esto implica que 
los caudales en el tiempo están tendiendo a 
mostrar valores más altos, y marca un cambio 
en los últimos 10 años, lo que podría asociarse 
con un aumento de las temperaturas. 

Es importante señalar que el método 
de análisis probabilístico utilizado para 
determinar el comportamiento de los caudales 
permite comparar, de acuerdo con criterios de 
significancia estadística, los valores máximos 
probables asociados con distintos periodos 
de retorno (30, 50 años, etcétera), y que son 
obtenidos mediante el uso de una FDP que 
presenta la mejor calidad de ajuste, en este caso, 
Gumbel. Así, es posible comparar los valores 
probabilísticos asociados con los periodos de 
datos dentro de una serie histórica, lo cual es 
una diferencia respecto de los test de tendencia, 
los que si bien indican cómo es la tendencia de 
cada uno de los periodos seleccionados dentro 
de la serie, no explican los valores asociados 
con dicho análisis. Sin embargo, el análisis 
probabilístico comparativo permite detectar 
incrementos O decrementos porcentuales de 
los montos anuales de los caudales máximos 
(Apaydin et al., 2006; Munka et al., 2007; 
Angulo-Martínez y Beguería, 2009; Lee y Heo, 
2011). 


Cálculo de la probabilidad de los rangos de 
cada serie según la función de Gumbel 


Con el fin de corroborar lo evidenciado en las 
gráficas, se calculó la probabilidad y los valores 
que adoptaría la variable de interés (caudales 
máximos), para el periodo de retorno de 30 años 
en cada una de las estaciones analizadas. Para 
esto se consideró la función de distribución 


de probabilidad de Gumbel. Ónoz y Bayazit 
(2001), en su estudio de caudales mínimos en 
16 regiones de Europa, determinaron en forma 
general que la función de distribución de 
Gumbel es la más recomendable en su uso. 

Se mantuvo la división de cada serie en los 
periodos 1965-1981 y 1982-2006, y los periodos 
1970-1981 y 1982-2006, para las estaciones, en 
donde se contempla menor cantidad de datos. 
Después se procedió a calcular la probabilidad 
y los valores que adoptaría la variable para el 
periodo de retorno nombrado (30 años), para 
conocer estadísticamente si se corrobora el 
aumento de los caudales punta en los últimos 
años. 

En el cuadro 6 se presentan los valores de 
caudales máximos asociados con el periodo de 
retorno de treinta años, comparados para los 
periodos 1965-1981 y 1982-2006. 

De acuerdo con los resultados presentados 
en el cuadro 6, en ocho de las nueve estacio- 
nes analizadas gráfica y estadísticamente, se 
observa que para una misma probabilidad se 
encuentran mayores valores de caudal máximo 
en el segundo periodo. 

Así, se puede decir que la gran mayoría 
de las estaciones aumenta sus caudales en su 
segundo periodo entre un 20 y 80% (diferencia 
porcentual), con respecto al que le precede 
(cuadro 6). 

Existen dos casos particulares, la estación 
Río Mapocho en Rinconada de Maipú para la 
Región Metropolitana, y la estación Río Claro 
en El Valle para la Región de O'Higgins. La 
estación Río Mapocho en Rinconada de Maipú 
hace una estimación de los valores de caudal, 
que supera el primer periodo en cuatro veces 
aproximadamente, lo que llama la atención, 
considerando que las demás estaciones no 
llegan al doble. 

La otra estación que se diferencia de las 
demás es El Valle en la Región de O'Higgins, 
estación que presentó valores que demuestran 
una disminución de los caudales entre el 
primer y segundo periodo considerado, 
aunque esta disminución es muy pequeña y 
casi despreciable. 
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Cuadro 6. Caudales máximos asociados con el periodo de retorno de 30 años para los periodos considerados 


para las distintas estaciones. 


Estación *Periodo de Periodo Periodo 
retorno 1965-1981 1982-2006 Diferencia % 
P(x<X) (m/s) (m/s) 
El Manzano (Región Metropolitana) 30 776 1178 52 
cena Nieves (Región de 30 333 385 16 
Ancoa (Región del Maule) 30 547 791 45 
San Alfonso (Región Metropolitana) 30 424 674 59 
Quilamuta (Región de O”Higgins) 30 342 613 79 
Nirivilo (Región del Maule) 30 223 290 30 
Rda. Maipú (Región Metropolitana) 30 179 933 421 
El Valle (Región de O'Higgins) 30 584 563 4 
Quiriquina (Región del Maule) 30 2289 2742 20 


* Probabilidad de excedencia. 


Según los resultados presentados en 
el cálculo de la variable en estudio, para 
el periodo de retorno de treinta años se 
corroboraría un aumento de los caudales 
punta en la última década. 

Con base en las interrogantes que surgen 
acerca de los posibles procesos de cambio 
climático que se estarían generando en 
los últimos años, estos resultados podrían 
ayudar a corroborar lo descrito por algunos 
autores en la revisión bibliográfica pre- 
sentada, en donde se ha señalado que existe 
un aumento de las temperaturas y, por 
consiguiente, una menor reserva de nieve en 
la cordillera, además de un posible retroceso 
glaciar (Vicuña et al., 2011). 


Test U Mann-Whitney 


Además, y para visualizar si las series 
provienen de la misma población o son eventos 
independientes, se realizó la prueba U Mann- 
Whitney a las seis estaciones graficadas 
anteriormente, y se sumaron las siguientes 
series comparativas: 1965-1985/1986-2006 y 
1970-1980/1990-2000. 

Esta prueba no paramétrica que no 
requiere de los supuestos de normalidad 
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e igualdad de varianzas, determina si dos 
muestras independientes provienen de la 
misma población (Montgomery y Runger, 
1996; Mendenhall y Sincich, 1997) (cuadro 7). 

De acuerdo con los resultados obtenidos 
en el cuadro 7, se puede observar que, en 
general, y entre un 80 y 90%, las series de 
datos no tienen diferencias significativas. 
Esto podría explicarse por la fuerte influencia 
glaciar en dichas estaciones, pertenecientes a 
la Región Metropolitana, lo que contrastaría 
con el resto de las estaciones que, en general, 
son eminentemente nivopluviales y podrían 
explicar el fenómeno expuesto. Sin embargo, 
cabe destacar el comportamiento de la estación 
Río Mapocho en Rinconada de Maipú, en 
donde todas sus series son significativamente 
diferentes. Esto podría deberse a que en 
la última década se ha dado un explosivo 
aumento urbano en la zona, sumado a 
situaciones meteorológicas, como lo es la 
concentración de las precipitaciones. Estas 
situaciones incidirían de forma directa en el 
comportamiento de los caudales máximos, 
además de que dicha estación se encuentra 
a la salida de la cuenca (Hawley y Bledsoe, 
2011). 
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Cuadro 7. Resultados test U Mann-Whitney, para tres estaciones de las regiones Metropolitana, O'Higgins y del Maule. 


Serie 
Región Estación 1965-1981/1982-2006 | 1970-1980/1990-2000 | 1965-1985/1986-2006 | 1975-1990/1991-2006 
Resultado 
Río Maipo en El Manzano Acepta H, Acepta H, Acepta H, Acepta H, 
Metropolitana | Río Maipo en San Alfonso Acepta H, Acepta H, Acepta H, Acepta H, 
Río M h Ri d 
A on A Ñ Acepta H Acepta H, Acepta H, Acepta H 
de Maipú t 3 . k 
Estero Alhué en Quilamuta Acepta H, Acepta H, Acepta H, Acepta H, 
O'Higgins AA pos AS Acepta H, Acepta H, Acepta H, Acepta H, 
Nieves y 2 Y E 
Río Claro en El Valle Acepta H, Acepta H, Acepta H, Acepta H, 
Río Ancoa en El Morro Acepta H, Acepta H, Acepta H, Acepta H, 
Maule Río Purapel en Nirivilo Acepta H, Acepta H, Acepta H, Acepta H, 
Río Longaví en Quiriquina Acepta H, Acepta H, Acepta H, Acepta H, 
Donde: 
H,; las medianas de ambas poblaciones son idénticas. 
H.; las medianas de ambas poblaciones son diferentes. 
Conclusiones Recibido: 22/05/12 
Aceptado: 03/08/2012 
Según el análisis realizado en el presente 
estudio, se concluye que los caudales presentan Referencias 


una tendencia a registrar mayores valores en 
las últimas dos décadas (1982-2006). 

El aumento de estos caudales varía desde 
un 100% hasta un 14%, llegando incluso a 
diferencias del 900%, al compararse las series 
1965-1981 y 1982-2006. 

Por otro lado, el aumento de los caudales 
máximos anuales podría confirmar, desde el 
punto de vista hidrológico, procesos de cambio 
climático, lo que ha generado un mayor aporte 
de caudales, debido a que no se observan otros 
fenómenos que afecten mayormente a los 
caudales. 

Por tanto, la situación descrita podría 
explicarse por el aumento de las temperaturas y 
se Observaría un cambio en el comportamiento 
nival y glaciar. Está situación corrobora lo 
descrito en los estudios de Carrasco et al. (2005 
y 2008), quienes señalan el incremento de la 
isoterma 0” en la zona central de Chile. 
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Abstract 


PIZARRO, R., BALOCCHI, E, VERA, M., AGUILERA, A., MORALES, C., VALDÉS, R., 
SANGUESA, C., VALLEJOS, C., FUENTES, R., ABARZA, A. € OLIVARES, C. Influence 
of climate change on peak flow behavior in central Chile. Water Technology and Sciences 
(in Spanish). Vol. IV, No. 2, April-June, 2013, pp. 5-19. 


This study analyzed the influence of climate change on peak flow behavior in the metropolitan 
region, Libertador Bernardo O'Higgins region and the area of Maule, central Chile. To this 
end, the Gumbel's PDE was applied, with which the best results for fitting peak flows were 
obtained, considering return periods of 10, 20, 30, 40, and 50 years. Thanks to the time-series 
segmentation of data, we concluded that the flow values tended to increase during the most 
recent periods. This is illustrated with mean values, trend charts and fitting of probability 
distribution functions. Finally, it was graphically and statistically shown that during recent 
years there has been an increase in maximum annual flows, which could be confirmed through 
hydrological climate change processes which have increasingly contributed to the flows 
through decreased snowpacks and glacial retreat from increased temperatures. 


Keywords: peak flows, climate change, probability density function, Mediterranean zone. 
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